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We analyze the thermal Casimir-Polder potential experienced by a ground-state molecule in a 
planar cavity and investigate the prospects for using such a set-up for molecular guiding. The 
resonant atom-held interaction associated with this non-equilibrium situation manifests itself in 
oscillating, standing-wave components of the potential. While the respective potential wells are 
normally too shallow to be useful, they may be amplified by a highly reflecting cavity whose width 
equals a half-integer multiple of a particular molecular transition frequency. We find that with an 
ideal choice of molecule and the use of an ultra-high reflectivity Bragg mirror cavity, it may be 
possible to boost the potential by up to two orders of magnitude. We analytically derive the scaling 
of the potential depth as a function of reflectivity and analyze how it varies with temperature and 
molecular properties. It is also shown how the potential depth decreases for standing waves with 
a larger number of nodes. Finally, we investigate the lifetime of the molecular ground state in a 
thermal environment and find that it is not greatly influenced by the cavity and remains in the 
order of several seconds. 

PACS numbers: 34.35. +a, 12.20.-m, 42.50.Ct, 42.50.Nn 



I. INTRODUCTION 

Casimir-Polder (CP) forces are a particular example of 
dispersion forces, which arise due to the fluctuations of 
the quantized electromagnetic field These forces oc- 
cur between polarizable atoms or molecules and metallic 
or dielectric bodies and can be intuitively understood as 
the dipole-dipole force that arises from spontaneous and 
mutually correlated polarization of the atom or molecule 
and the matter comprising the body. Casimir-Polder 
forces at thermal equilibrium have been commonly in- 
vestigated in the linear-response formalism 0-|U • Studies 
of a wide range of geometries such as semi-infinite half 
spaces , thin plates [a @] , planar cavities H , spheres 
and cylinders [§] as well as cylindrical shells [a IS US] have 
revealed that thermal CP forces are typically attractive 
in the absence of magnetic effects. 

Recent theoretical predictions filfl as well as experi- 
mental realizations [12|] for CP forces in thermal non- 
equilibrium situations have pointed towards interesting 
effects which arise when an atom at equilibrium with its 
local environment interacts with a body held at different 
temperature. In particular, depending on the tempera- 
tures of the macroscopic body and the environment, the 
force can change its character from being attractive to 
repulsive and vice versa. 

Non-equilibrium between atom and local environment 
can be investigated by means of normal-mode quantum 
electrodynamics (QED) [lU, [l4j] or macroscopic QED in 
absorbing and dispersing media [l5|, [IS]- In this case, 
thermal excitation and de-excitation processes lead to 
resonant contributions to the force. As discussed in a 



recent Letter by two of the present authors [17[ (cf. sim- 
ilar findings reported in Ref. [3]), these resonant forces 
produce a different force from that obtained through the 
standard approach, a perturbative expansion of Lifshitz' 
formula using the ground state polarizability of the atom. 
Only when the atom is fully thermalized, i.e., when it is 
in a superposition of energy states as given by the Boltz- 
mann distribution, do the two approaches yield the same 
result, and only when the correct, thermal polarizability 
is employed. For most atoms the resonant contribution is 
small because the respective excitation energies are much 
larger than the thermal energies, hence the atom is es- 
sentially always in its ground state. 

For diatomic polar molecules such as LiH or YbF, how- 
ever, whose lowest rovibrational eigenstates are typically 
separated by energies that are small on a thermal scale, 
the situation is changed, and the thermalized CP force 
can differ drastically from the standard 'Lifshitz-like' ex- 
pression. An investigation into these effects was under- 
taken in Rcf. where it was found that for YbF out- 
side a metallic half-space the fully thermalized CP force 
is smaller than the non-resonant force alone by a fac- 
tor of 870. These results could be of importance for the 
trapping of Stark-decelerated polar molecules near 
macroscopic bodies. 

Equally interesting is the observation that even for an 
atom or a molecule in its ground state the resonant part 
of the Casimir-Polder force has a long-range and spatially 
oscillating contribution, due to propagating modes [1J, 
19]. While this oscillatory behavior dies out as the system 
thermalizes, the thermalization time of a ground state 
molecule can be quite long, often several seconds [2l| . 
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The oscillating propagating potential reported in [19(, 
unfortunately, was found to be too small in amplitude 
to be useful for guiding purposes, but nonetheless points 
to interesting applications if a way could be found to 
enhance these oscillations. 

In the present article we investigate the use of a pla- 
nar cavity to enhance the amplitude of the potential os- 
cillations. This geometry has been discussed in detail in 
conjunction with excited atoms in a zero temperature en- 
vironment, where an oscillating, standing- wave potential 
is known to occur [22T - I281 ] . For ground-state molecules 
in a cavity at finite temperature, an enhancement of up 
to two orders of magnitude will indeed be shown to be 
possible when the cavity width is fixed at the resonant 
length a = ttHc/Eiq where Eiq is the energy separation of 
the ground state and first excited state of the molecule. 

The paper is organized as follows: The general formal- 
ism of the CP force on a molecule in a cavity is developed 
in Sec.|Tll and numerical calculations for a gold cavity are 
undertaken in Sec. Ill Al where we also show how the po- 
tential is enhanced as the cavity approaches the resonant 
width. In Sees. Ill Bl and III CI we investigate strategies 
for further enhancing the potential by considering differ- 
ent cavity resonances and molecular species. The scaling 
of the potential with the reflectivity of the cavity is in- 
vestigated numerically and analytically in Sec. Ill Dt and 
we thereafter discuss how further enhancement can be 
achieved using a cavity of parallel Bragg mirrors tuned to 
frequency uiqi = E i/h and normal incidence fSec. Ill El) . 
Finally, in Sec. Ill Fl we investigate the effect of the cavity 
on the thermalization time of a molecule initially pre- 
pared in its ground state and find that this remains in 
the same order of magnitude as in free space, typically in 
the order of seconds. We summarize our result in Sec. lIIII 
and provide a guide to further investigations. 



II. THERMAL CASIMIR-POLDER POTENTIAL 
IN A PLANAR CAVITY 



We consider a polar molecule with energy eigenstates 
|n), eigenenergies hio n , transition frequencies uj mn — 
U! m — Lu n and dipole matrix elements d m „, which is pre- 
pared in an incoherent superposition of its energy eigen- 
states with probabilities p n . As shown in Ref. [17[, 
the CP force is conservative in the perturbative limit, 
F(r) = — V£7(r), where the associated CP potential is 
given by 



U(r) = -J2pnU n (r), 



(1) 



and the potential components for an isotropic molecule 
read 

oc 

U n (r) = f i k B TY / '^^ n (iQReTr G«(r, r, i£) 
+ y y^^ fc |d» fc | 2 {8(a; fc „)n(a; fc „) 

k 

- Q(u> nk )[n(u nk ) + l]}ReTrG (1) (r,r, |w nfc |). 

(2) 

where (1q is the free-space permeability, fee is Boltz- 
mann's constant, £j = 27tj7cbT/7i is the jth Matsubara 
frequency, and (r, r', w) is the scattering part of the 
classical Green tensor of the geometry the molecule is 
placed in. The prime on the Matsubara sum indicates 
that the j = term is to be taken with half weight. The 
molecular polarizability is given by 



a„(w) = lim 



Ink 



d„fe 



w - Lu kn + ie 



(3) 



The photon number follows the Bose-Einstein distribu- 
tion, 



exp 



- 1 



(4) 



The first sum in Eq. ^) is the non-resonant force, rem- 
iniscent of that obtained by a dilute-gas expansion of 
Lifshitz' formula Q. The second sum is the resonant 
contribution to the force. We will see how it splits natu- 
rally into a propagating plus an evanescent part. 

We assume the molecule to be placed within an empty 
planar cavity bounded by two identical plates of infinite 
lateral extension with plane parallel surfaces, separated 
by a distance a. We choose the coordinate system such 
that the cavity walls are normal to the z axis at z = ±a/2 
(z = being the center of the cavity) and denote direc- 
tions in the xy plane by the symbol _L. The scatter- 
ing Green tensor of the system is well known (cf., e.g., 
Ref. (29|), and the relevant diagonal elements inside the 
cavity are given by 



G^J(z,z;u,kx) 



^^- e ^ C082 f3z, 
z lPa cos2/3z, 



P D 
uj 2 [3 D p 



cos 2f3z, 



(5a) 
(5b) 
(5c) 



•(r-r')j 



where we have performed a Weyl expansion, 

G (1) (r,r»= / ^GW( Z , Z ',k llW ) e * 

(6) 

taken the coincidence limit r — » r' and dropped all 
position-independent terms (which give rise to an irrel- 
evant constant contribution to the CP potential). Here, 
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r s , r p are the reflection coefficients of the (identical) cav- 
ity walls for s,p polarized waves and we have defined 



D„ = l- Ke 



2 2i/3a 



(3 = yjw 2 /^ -k\ 



(7) 

(8) 



The square root is to be taken such that lm/3 > 0. 
When the cavity walls are homogeneous, semi-infinite 
half-spaces of an electric material of permittivity e(w), 
the reflection coefficients can be written simply as 



p- 


v7? 2 + 


-(e- l)w 2 /c 2 


+ 




(e- l)w 2 /c 2 


e/3 




+■ (e- l)w 2 /c 2 


ep 




+■ (e- l)w 2 /c 2 



(9a) 
(9b) 



where again the square roots are chosen such that their 
imaginary part is positive. 

Adding Eqs. (|5a |) ~ ([5c|) and partially performing the 
Fourier integral by introducing polar coordinates in the 
xy plane, the trace of the Green tensor of the cavity reads 



TrG 



i 

27d 



(r,r,w) 

k i dk i 



/3 



c 2 /3 2 rp_ 

LO 2 D„ 



E 



D n 



i{3a 



co$2/3z. 
(10) 



This result can be substituted into Eq. (j2J to obtain the 
thermal CP potential of a molecule in an arbitrary in- 
coherent internal state. In the following, we will assume 
the molecule to be prepared in its ground state, so that 
the thermal CP potential is given by 

oo 

U(r) = » k B T '$a(i£j)Be Tr G (1) (r, r, »&) 

3=0 

+ y E < n M \d k | 2 Re Tr G« (r, r, Uk0 ), 



(11) 



[a(uS) = ao(w), ground-state polarizability] together with 
Eq. ((TU)) . The first term is the non-resonant part of the 
potential, it depends on the Green tensor taken at purely 
imaginary frequencies. Since ft is purely imaginary in this 
case, the Green tensor (fTUj) and hence the non-resonant 
potential is non-oscillating as a function of position. The 
second term in the CP potential is the resonant contri- 
bution, which depends on the Green tensor taken at real 
frequencies. The integral over k± in this case naturally 
splits into a region < k± < uo n k of propagating waves 
in which /3 is real and positive, and a region ui n k < fcj_ 
of evanescent waves in which f3 is purely imaginary. The 
contributions from propagating waves are oscillating as a 
function of position due to the term cos 2(3 z, while those 
from evanescent waves are non-oscillating, just like the 
non-resonant part of the potential. The total poten- 
tial (jlip can thus be separated into non-resonant (first 
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FIG. 1: Casimir-Polder potential of a ground-state LiH 
molecule inside a gold cavity of width a — 500/im at room 
temperature (T = 300K). The non-resonant, propagating, 
and evanescent contributions to the total potential are shown 
separately. 



term), propagating (contributions to second term with 
< k±_ < uj n k), and evanescent components (contribu- 
tions to second term with uj„k < k±) according to 



U(z) = U nl (z) + U pr (z) + U cv {z) 



(12) 



To illustrate the behavior of the total potential and 
its three components, we consider a LiH molecule in its 
electronic and rovibrational ground state placed inside a 
gold cavity. The permittivity of the (semi-infinite) cavity 
walls may be computed using the Drude model 



i 



uj(u + 17) 



(13) 



with w p = 1.37 x 10 16 rad/s and 7 = 5.32 x 10 13 rad/s 
[30l | . As shown in Ref . [l!| , the CP potential of ground- 
state LiH is dominated by contributions from the ro- 
tational transitions to the first excited manifold, with 
the respective transition frequency and dipole matrix el- 
ements being given by uj k = 2.78973 x 10 12 rad/s and 
J2 k |d 0fc | 2 = 3.847 x lCr 58 C 2 m 2 , respectively E'l . The 
potential (fTTj) and its three components (fT2"j) for a cavity 
of length a — 500/im at room temperature (T = 300K) is 
shown in Fig. Q] as the result of a numerical integration, 
where Eqs. ©-([TU]) have been used. For transparency, 
we have shifted all three components such that they van- 
ish at the center of the cavity. It is seen that the non- 
resonant potential is attractive and has a maximum at 
the center of the cavity, while the evanescent potential 
is repulsive and has a minimum at the cavity center. As 
in the case of a single surface [l^| these two contribu- 
tions partially cancel, where the attractive non-resonant 
contribution is slightly larger and leads to an attractive 
total potential in the vicinity of the cavity walls. The 
propagating part of the potential is spatially oscillating 
and finite at the cavity walls, it dominates in the central 
region of the cavity where it gives rise to well-pronounced 
maxima and minima. 
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It is natural to wonder whether these potential min- 
ima might be used for the purpose of guiding of polar 
molecules. With this in mind, we will in the following 
discuss strategies of enhancing the depth of the potential 
well by analyzing the dependence of the potential on the 
molecular species as well as the geometric and material 
parameters of the cavity. 



A. Cavity-induced enhancement of the potential 

We begin our analysis by discussing the dependence of 
the potential on the cavity width. The one-dimensional 
confinement of the propagating modes in a cavity of 
highly reflecting mirrors leads to the formation of stand- 
ing waves and associated cavity resonances. When the 
molecular transition frequency coincides with one of these 
resonances, the thermal CP potential can be strongly 
enhanced: When the squared reflection coefficient is 
close to unity, the denominator D a of Eq. ([7|), featur- 
ing in the Green tensor, becomes small if the exponential 
exp(2i/3a) is equal to unity, resulting in a strong enhance- 
ment of the potential U-p X . This happens for normal inci- 
dence (fc_L = 0) of the propagating waves, when the res- 
onance condition 2u>Qka/c = 2iri>, v G N is fulfilled. In 
other words, the cavity length has to be equal to a half- 
integer multiple of the molecular transition wavelength 
A/co = 2nc/uj k0 : 

a = uX k0 /2, i/£N. (14) 

We say that the molecular transition coincides with the 
vth cavity resonance. 

The cavity-induced enhancement of the thermal CP 
potential is illustrated in Fig. [31 where we show the total 
thermal CP potential of a ground-state LiH molecule in 
gold cavities of widths such that the molecular transition 
is close to the second cavity resonance Xko- As seen, the 
amplitude of the spatial oscillations, associated with the 
propagating part of the potential U pr , sharply increases 
as the cavity width approaches Xko- For comparison, 
we have also displayed the potential of a single plate at 
z = — Afeo/2, where Eq. (TTU)) for the cavity Green tensor 
has been replaced with the single-plate result jjjl 

TrG (1) (r,r',w) 

_ j_ r°° k ± dk ± I" _ 2 <?fP ] ip( a +2z) 

U (7 — S,p L J 

The comparison shows that the amplitude of the oscilla- 
tions, while hardly visible for the single plate, is strongly 
enhanced for a cavity. The depth of the potential mini- 
mum at the center of the cavity with respect to the neigh- 
boring maxima is increased by a factor 6.7 when using a 
resonant cavity rather than a single plate. 
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FIG. 2: Cavity-induced enhancement of the thermal Casimir- 
Polder potential of a ground-state LiH molecule inside gold 
cavities of various widths close to the second resonance a = 
Afco = 673/xm. The potential of a single plate at z — — Afco/2 
is also displayed. 
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FIG. 3: Propagating part of the thermal Casimir-Polder po- 
tentials of a ground-state LiH molecule inside gold cavities of 
widths a = v\ko/2 corresponding to resonances of different 
orders (y = 1, 2, 3). 



B. Different cavity resonances 

In the following, we are interested in the cavity- 
enhanced oscillations of the thermal potential. As seen 
from Fig. [21 they set in at some distance away from the 
cavity walls where the potential is well approximated by 
its propagating- wave contribution U pl . We can therefore 
restrict our attention to this part of the total CP poten- 
tial. The (propagating-wave) potentials associated with 
different cavity resonances v are shown in Fig. [31 It is 
seen that the order v of the resonance corresponds to the 
number of maxima of the potential. Potentials associated 
with resonances of order v > 2 have minima. The am- 
plitudes of the oscillations become generally smaller for 
higher resonance orders v. As seen from the case v = 3, 



FIG. 4: Depth AU of the v = 2 potential minimum for ro- 
tational and vibrational transitions of various ground-state 
polar molecules inside gold cavities at T = 300K. The differ- 
ent non-degenerate transitions of OH and OD are labeled as 
(a)-(d) in order of ascending frequencies, cf. Ref. pl| . 



Frequency (rad/s) 

FIG. 5: Frequency-dependence of the depth AU of the v = 2 
potential minimum for polar molecules inside gold cavities at 
T = 300K. 



the minima and maxima are slightly more pronounced 
towards the cavity walls. 

The scaling of the potential minima with the resonance 
order as observed in Fig. [3] can be confirmed by an analyt- 
ical analysis. For each cavity order v, we define AU V to be 
the depth of the deepest potential minimum with respect 
to the neighboring maxima. As suggested by Fig. [3J this 
deepest minimum will always be the one closest to the 
cavity walls. Cavity QED problems can often be solved 
analytically under the simplifying assumption that reflec- 
tion coefficients are independent of the transverse wave 
number k± [3l], [12], and this method is also successful 
here. As shown in App.|3 in the perfect conductor limit 
r p = —r s = r — > 1, we have the simple scaling law 



AU V 



1 



(16) 



For imperfect conductors, the AU V will decrease some- 
what less slowly with v. 

The analytical scaling law obtained on the basis of sim- 
plifying assumptions supports the observation from the 
numerical results in Fig. [3J that the v = 2 resonance pro- 
vides the deepest potential minimum. In view of poten- 
tial guiding, we can therefore restrict our attention to 
this case, AU = A£/ 2 . 



C. Different molecular species 

The CP potential depends on the molecular transition 
in question via the respective transition frequencies and 
dipole matrix elements. Using the molecular data as 
listed in Ref. [2l[, we have calculated the depth of the 
v = 2 potential minimum for both rotational and vibra- 
tional transitions of the polar molecules LiH, NH, OH, 
OD, CaF, BaF, YbF, LiRb, NaRb, KRb, LiCs, NaCs, 
KCs, and RbCs; the results are displayed in order of de- 
scending AU in Fig. [4] The figure shows that the deepest 
potential minima are realized when using the rotational 
transition of LiH, the vibrational transition of BaF or the 
dominant rotational transition of OH, followed by YbF 



(vibrational), LiRb (vibrational), NH (rotational), OD 
(dominant rotational transition), and CaF (vibrational). 

The variation in the depth for different molecules is 
partly due to its dependence on the molecular transi- 
tion frequency. As shown in Ref. (l9| . the resonant 
part of the CP potential of a single plate is propor- 
tional to uj^n^ko) for a good conductor with frequency- 
independent reflectivities. This remains true in the case 
of a cavity. In addition, the amplitude of the oscilla- 
tions is inversely proportional to the molecule-wall sepa- 
ration. The largest potential maximum being situated at 
z — a/2 = oc 1/uiko, its height carries an additional 

Wfco-proportionality. The dependence of the potential- 
minimum depth on molecular transition frequency can 
thus be given as 

att 3 / s j w feo for ftwfco < fc B T, 
AU V cx uJi n n(uJko) oc < V . 

fc0 V ; | e -ft<WteT for /i Wfc o > fc B T. 

(17) 

As shown in Sec. Ill D[ this scaling becomes exact for 
cavities with frequency- and ^-independent reflectivi- 
ties. For real conductors, the decrease of AU V for high 
frequencies will be stronger than given in Eq. (JTTJ) due 
to the decrease of the reflection coefficients. Note that 
Eq. (fT7|) also shows that AU V becomes larger for higher 
temperatures due to the increased thermal-photon num- 
ber. Again, this only holds when disregarding the tem- 
perature dependence of the reflection coefficients, cf. also 
Sec. |HE] below. 

The frequency-dependence of AU V is illustrated in 
Fig. [5] where we have plotted its values normalized by 
dividing by the transition dipole moments d 2 (d 2 = 
Ylk |dofc| 2 )- The transition frequencies of some of the 
molecules investigated are indicated in the figure. In 
particular, the vibrational transitions of BaF and YbF, 
which have been seen in Fig. 2] to give rise to large 
potential-minimum depths, are very close to the peak 
of the function u)%. Q n (cj ^o) , which is at cuko = 1-H X 
10 rad/s for room temperature. 

The other main dependence of AU V on the molecular 



G 



species and transition is the proportionality to the mod- 
ulus squared of the transition-dipolc moments, 



Ok- 1 



The transition-dipole moments are typically larger for 
rotational transitions than for vibrational ones. For this 
reason, the rotational transition of LiH gives rise to the 
largest minimum depth although the vibrational transi- 
tion frequencies of BaF and YbF are much closer to the 
peak frequency 1.11 x 10 14 rad/s. 



D. Scaling with reflectivity 

The cavity-induced enhancement of the thermal CP 
force strongly depends on the reflectivity of the cavity 
walls. To understand this dependence in more detail, let 
us for simplicity investigate how the height of the single 
maximum for a v — 1 resonance depends on reflectivity. 
The scaling of the potential extrema with reflectivity is 
the same for all v as is shown in App. 13 so considering 
the simplest case will suffice. 

We begin by writing the propagating part of the reso- 
nant CP potential associated with a single transition in 
the form 

U pi (z) = ~n(u ko )\d 01 \ 2 I(<fi) (19) 

6Eq 



where we have introduced the dimcnsionless position 



and the integral 



!(</)) = Im 



,ko/c dkj_k± 
2tt/3 



z 
a 



2(3' 



(20) 



<4o 



x e ipa cos2/3a<t). (21) 



As in Sec. Ill B[ we consider the simple model case 
of frequency- and fcj_ -independent reflection coefficients 
r„ = — r s = r. With this assumption, 



0. 



(22) 



After introducing the dimensionless integration variable 
x = 2/3a with k±dkj_ = —Aa 2 xdx, the integral above 
takes the form 



87ra 3 



Im 



dx 



p 2 j-jtcc / 2 



COS 0X 



1 



(23) 



where xq = 2ukoa/c. For the v = 1 resonance, we have 
a = Aofc/2 = nc/uj k o, so x = 2tt. 

The required height of the potential maximum at the 
cavity center (z = 0) with respect to the value of U pT at 



the cavity walls z = ±a/2 is proportional to the difference 
1(0) We have 



(18) /( i/2) = 



2tt 

Im / dxx 2 
lGna 6 .In 1 — r z e tx 



167T0 3 



Im 



dxx 2 [l + (l + r- 2 )Li (r 2 e M )] 

(24) 



where the polylogarithmic function is defined as 

oc u 



fc=i 



(25) 



The first term in Eq. ((2"4"]) is real and does not contribute. 
The second one is easily calculated using the relation 



JdzU s (Ae bz ) = ^U s+1 (At 



C 



(26) 



valid for arbitrary constants A, b where \A\ < 1. Partially 
integrating this relation twice and substituting the result 
for A — r 2 , b = i into Eq. (f24| , one finds 

r -A- r~ 1 f 4-tt 2 
1(1/2) =-?—l Ya \—U 1 (r 2 )+ATrU 2 (r 2 ) 



(27) 



167ra 3 i 
7r(r + r _1 ) 



4a 3 



ln(l 



where we have noted that hii(z) = — ln(l — z). In the 
limit of high reflectivity, S = 1 — r — > 0+ this exact result 
shows the asymptotic behavior 

1(1/2) - ^ (In 5 + In 2) for 6 -> 0+, (28) 

with the first correction term being of order S. 

The calculation of 1(0) is only slightly more involved. 
We have 

/(0)=„lm/ dx 



87ra 3 



87ra 3 



1 — r 2 e ix 

2k 



^r 2fe Im / dxx 
1=0 Jo 



(29) 



By partial integration we obtain 



Im 



2tt 



4tt 2 



(l + \) (l+\f 



(30) 



dxx 2 e lx ^ 

10 \- i 2 ' 

After substitution of this result, the sum over I can be 
performed by using the relations (cf. §1.513 in Ref. [33| ) 

r-22 I ]__!_»• 



£ 



In 
r 1 



ln<5 + ln2 for 5^0+ (31) 



(leading corrections being of the order 5\nS) and 



E 



1=0 



X.) j X j 

E/3 "E (27)3 



1=1 



7C(3) for 6^0+, 
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where £(z) is the Riemann zeta function. We thus find 

J(0)~-/ifM-ln2+^C(3)] for 5^0+, (32) 

with the first correction again being of order 5 In 5. 

Substituting the results and (J35J into Eq. (O, the 
difference between the maximum and minimum values of 
the v = 1 propagating potential reads 



M0)-EV(a/2) 
for 5 -> 0+. 



Tr^fc |d ofc | 2 rt(a; fc o) 



ln5- 



7C(3) 
2tt 2 

(33) 



This result being representative of the case of arbitrary 
v, we can conclude that 



AU U cx ln(l - r) 



(34) 



in the limit r — > 1. In the case where the reflection coef- 
ficients are not the same for both polarizations but still 
assumed constant, the coefficient of the term oc ln<5 in 
Eq. ([33)1 will change, leading to a slight quantitative but 
no qualitative difference to the scaling of the potential- 
minimum depth. Note that Eq. (|33|l immediately implies 
the scaling law (fT7|) for the frequency dependence of AE/„. 

The fact that the potential depth diverges only log- 
arithmically as reflectivity tends to unity poses severe 
restrictions on the potential which is obtainable using 
a planar cavity. The mathematical reason for the rela- 
tive weakness of the resonance is that the integrand of 
the fcj_-integral only becomes large at a single point, at 
k± = 0. The physical reason is that the photonic modes 
in the cavity are only confined in one out of three spa- 
tial dimensions. We conjecture that the potential due to 
the resonant CP force on a ground state molecule can be 
much increased by a resonant cavity if confinement is im- 
posed in two or even three dimensions, i.e. in a cylindrical 
or spherical cavity. 

The logarithmic scaling law of the potential-minimum 
depth Af7„ for the case v — 2 is confirmed by a nu- 
merical calculation in which reflection coefficients are set 
constant, r p = —r s = r and close to unity. The result 
for the rotational transition of LiH is shown in Fig. 
where the exact result for a gold cavity is also included. 
By comparing the latter curve to the potentials for con- 
stant reflection coefficients, one can read off the relatively 
small 'effective' reflectivity of gold between I — I0~ 2 and 
1 — 10~ 3 at the respective transition frequency of LiH. 
For a molecule with a smaller eigenfrequency uiko the 
gold cavity does slightly better because the permittivity 
is larger. Consider the vibrational transition of YbF with 
Wko ~ 9 • 10 10 rad/s as an example, for which the 'effec- 
tive' reflectivity of the gold cavity (in the sense of figure 
in]) increases to about 1 — 10~ 3 5 . 



E. Enhanced reflectivity using Bragg mirrors 

In contrast with the non-resonant CP force which de- 
pends on a very broad band of frequencies, the resonant 
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FIG. 6: Propagating part of the thermal CP potential for a 
ground-state LiH molecule inside different cavities with con- 
stant reflection coefficients. The rotational transition of LiH 
is assumed to coincide with the v = 2 resonance of the cavi- 
ties. For comparison, the exact result for a gold cavity is also 
shown. 



part of the ground state force on a two-level molecule 
depends on the reflection properties of the cavity at a 
single frequency, lj — Uko- In addition, the resonance of 
the cavity is also associated with a single value of the 
wave vector k^, namely normal incidence. An enhance- 
ment of the propagating potential hence does not require 
a good conductor like gold which is a good reflector for 
a broad range of frequencies and all angles of incidence; 
instead, cavity walls whose reflectivity has a sharp peak 
at normal incidence and the single frequency uJko are suf- 
ficient. The obvious candidate is to use multilayer Bragg 
mirrors, which consist of alternating layers of two differ- 
ent materials, each layer of thickness being equal to one 
quarter of the wavelength Aio = 2i:/nujkQ in that layer 
where n is the respective refractive index. 

The reflection coefficient of a stack of layers with per- 
mittivities Sj and thicknesses dj is found by recursive use 
of the formula 



Ti jk- 



Uj +rj k (i..-) 



1 + r ij r jfc(/--) e2l/3jd3 



(35) 



which relates the reflection co- 



(ft = ^/n^/c 2 

efficient of a set of three adjacent layers ijk ■ ■ ■ (and all 
the layers behind) to the respective result for the next set 
of adjacent layers jkl ■ ■ ■ . If the fcth layer is the last one 
of the stack, the coefficients rj^n...) reduce to the two- 
layer coefficients rjk- In straightforward generalization 
of Eqs. ([9a| and (|9b|) . the two-layer coefficients read 



A - ft 



ft + ft 

£oft + £ift 



(36a) 
(36b) 



for s- and p-polarized waves, respectively. The Casimir 
effect for such multilayer stacks has been extensively 
studied in the past [13-37]. 



-Dissipation as indicated 

-No dissipation 

- 3m{e} reduced by factor It 




FIG. 7: Reflection coefficients of Bragg mirrors for normal 
incidence at the rotational transition frequency of LiH vs the 
number of double layers JV. Above: GaAs/AlAs mirror, where 
the results for reduced and vanishing absorption are also dis- 
played for comparison. Below: Vacuum/sapphire mirror at 
two different temperatures T — 77K and 300K. 



A very common pair of materials to use for Bragg mir- 
rors is GaAs and AlAs. At the rotational transition fre- 
quency of LiH, the permittivity of the two materials can 
be roughly given as eqaAs = 12.96 + 0.02i [H, HH and 
saiAs = 10.96 + 0.02i [40f. The reflection coefficient of 
a GaAs/ AlAs Bragg mirror is plotted as a function of 
the number of (double) layers TV in the upper panel of 
Fig. [3 For a given N, the Bragg mirror consists of 2AT+1 
layers in total, i.e. N pairs of GaAs and AlAs layers of 
thickness Aio/4 (beginning with GaAs) and a terminat- 
ing GaAs layer of infinite thickness. As Fig. [7] shows, 
the reflectivity initially increases for increasing TV and 
then eventually saturates for N > 30 to some finite value 
where 1 — Rer ~ 10 -2 . This saturation is due to ab- 
sorption, as is illustrated by the other two curves, where 
we have given the results that would be obtained for a 
reduced or vanishing imaginary part of the permittivi- 
ties. For a reduced imaginary part, the saturation sets 
in for higher N, and consequently to a lower 6. In the 
absence of absorption, the reflectivity could be brought 
arbitrarily close to unity by adding more and more layers. 

A higher reflectivity could hence be obtained by using 
materials with very small dielectric loss. One example of 
such a Bragg mirror could be alternating layers of vacuum 
and sapphire, which can have an extremely low loss tan- 
gent (Ime/Ree ~ 10~ 5 and 10~ 7 at room temperature 
and 77K, respectively [4l|) combined with a refractive 
index considerably larger than unity (Ree ~ 10 [421]). 
Using the approximative values £ sap ph = 10 + 10~ 4 i at 
300K and £ sa pph = 10 + 10~ 6 i at 77K, we have computed 




Position z (jam) 

FIG. 8: Resonant part of the thermal CP potential associated 
with the rotational transitions of a ground-state LiH molecule 
at v — 2 resonance with a gold cavity and a cavity bounded 
by vacuum/sapphire Bragg mirrors at two different tempera- 
tures: 77K above, 300K below. The solid black lines represent 
calculations at constant reflection coefficients as in Fig. [6] the 
corresponding values of 1 — r decrease in powers of 10 from 
10 -2 (lowest curve) to 10 -7 (highest curve). The same per- 
mittivity is used for gold for both temperatures. 



the reflection coefficients of the vacuum/sapphire mirror 
as displayed in the lower panel of Fig. [7] At room temper- 
ature, the coefficient saturates atiV>6to6 = 5.5T0~ 6 . 
At T = 77K, the reflection coefficient saturates at N > 8 
to S = 5.5 ■ 10~ 8 , the the increase in reflectivity is ob- 
viously due to the reduction of material absorption for 
the lower temperature. Note that in comparison to the 
GaAs/ AlAs mirror, the number of layers required for sat- 
uration is significantly lower because of the larger dielec- 
tric contrast; and the room-temperature reflectivity at 
saturation is increased by about four orders of magni- 
tude. 

The resulting propagating part of the resonant CP po- 
tential at resonant cavity width using the sapphire/vac- 
uum Bragg mirror at T = 77K and 300K are shown in 
Fig. [SI where the corresponding graphs at various con- 
stant reflection coefficients have also been displayed for 
reference. The 'effective' reflection coefficients achieved 
at the two temperatures are around 6 = 10~ 4 ' 8 and 
S = 10~ 6,7 respectively, and the potential depths ap- 
proximately a factor 2.45 and 1.77 greater than that of 
the gold cavity at the same temperatures. Note, how- 
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ever, that the effect of the enhanced reflectivity at 77K 
is counteracted by the overall decrease of the potential 
due to the lower photon number. 



F. Lifetime of the ground state in the cavity 

Resonant CP potentials are only present for molecules 
which are not at equilibrium with their thermal environ- 
ment, i.e., on a time scale given by the inverse heating 
rate [19|] . When enhancing the thermal CP potential via 
a resonant cavity, it is necessary to ascertain that the si- 
multaneous cavity-enhancement of heating rates does not 
reduce the lifetime of the resonant potential by so much 
as to render it experimentally inaccessible. We show in 
the following that the lifetime of the molecular ground 
state is not radically changed even by the presence of a 
resonant planar cavity. 

The total heating rate of an isotropic molecule out of 
its ground state may be written as [17;] T — Tq + r cav 
where 



r„ 



EfcMofcl 2 w% n(w k o) 
3irhc 3 Eo 



(37) 



is the heating rate in free space and 



2^o 
3h 



E 

k 



d fc| 2 Wfc ', 



i(cj fc0 )ImTrG (1) (r,r, 



is its change due to the presence of the cavity. Apart 
from the prefactor, this additional term has the same 
form as the expression for the potential, except that the 
imaginary part of the Green tensor is taken rather than 
the real part. 

In Sec. Ill Dl we had shown that for real and constant 
reflection coefficients, the Green tensor exhibits a loga- 
rithmic divergence as r — > 1 with a purely real coefficient, 
whereas all other contributions remain finite. This shows 
that the imaginary part of the Green tensor responsible 
for the decay rate can be expected to remain finite even 
for strongly increased reflectivity. It follows that the pres- 
ence of the cavity does not drastically change the lifetime 
of the ground state of the molecule, which will typically 
be in the order of seconds. This is confirmed in Fig. [S] 
where we display the ground-state heating rate of a LiH 
molecule inside a v — 1 gold cavity and near a gold half- 
space. The lifetime is reduced by only a factor 2 at the 
center of the cavity, remaining in the order of seconds. 



III. CONCLUSIONS AND OUTLOOK 

We have studied the thermal Casimir-Polder potential 
of ground-state polar molecules placed within a planar 
cavity at room temperature. As was previously found 
in Ref. [HI, the resonant absorption of thermal photons 
by a molecule gives rise to spatial oscillations of that 
potential. Our results demonstrate that the amplitude of 
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FIG. 9: Ground-state rotational heating rate of a LiH 
molecule in a gold cavity at v = 1 resonance a — -Kc/uJho 
(solid line) and to the right of a gold half-space (dashed line) 
at room temperature (T = 300K). 



these oscillations is enhanced when placing the molecule 
inside a suitable cavity such that a molecular transition 
frequency coincides with a cavity resonance. We have 
analyzed the dependence of this oscillating potential on 
the parameters of the molecule and the cavity by both 
analytical and numerical means and found that the depth 
of potential minima . . . 

• Cavity resonance: . . . decreases with increasing or- 
der of the cavity resonance approximately as l/v; 



Molecular eigenfrequency: 



is proportional to 



u>k n(uika) for good conductors, where n(uJko) is the 
thermal photon number; 

• Molecular dipole moments: ... is proportional to 
the modulus squared ^ fc |do/c| 2 = d 2 of the total 
transition dipole moment; 

• Temperature: . . . increases with temperature due to 
an increase of the thermal photon number n(wfco); 



• Reflectivity of cavity walls: 
for high reflectivity r. 



scales as ln(l — r) 



In view of observing this potential and possibly uti- 
lizing it for the guiding of cold polar molecules, these 
observations imply the following strategies for enhancing 
the depth of the potential minima: 

• Cavity resonance: The v = 2 resonance is most 
suitable, since it gives the deepest minimum. 

• Molecular species: At room temperature, the deep- 
est minima are realized for molecules whose tran- 
sitions are not too far from the peak frequency 
1.11 x 10 14 rad/s and which at the same time feature 
suitably large transition dipole moments. Good 
candidates are, e.g., LiH (rotational transitions), 
BaF (vibrational transitions) or OH (rotational 
transitions). 
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• Cavity walls: Highly reflecting cavities are required 
in order to enhance the potential. Bragg mirrors 
consisting of materials with small absorption such 
as sapphire are favorable to single layers of good 
conductors like gold. 

• Temperature: Temperatures should be in the range 
of room temperature or even higher in order to 
achieve large photon numbers. This should be bal- 
anced, however, against the adverse reduction of 
reflectivity of most materials with increasing tem- 
perature. 

With an optimum choice of all these parameters, the 
planar cavity can be used to enhance the resonant po- 
tential by one or at most two orders of magnitude with 
respect to the single-plate case. However, the thermal 
potentials achievable with planar cavities are in all like- 
lihood still too small to facilitate the guiding of polar 
molecules. 

The limitations of the enhancement of the potential 
in a planar cavity are ultimately due to the weak (loga- 
rithmic) scaling with reflectivity. A stronger scaling may 
be expected in geometries providing mode confinement 
in more than just one dimension such as cylindrical or 
spherical cavities. This will be investigated in a future 
publication. Note that apart from the different expected 
scaling with reflectivity, all other conclusions regarding 
the dependence of the potential on the relevant molecular 
and material parameters as given above hold irrespective 
of the geometry under consideration. The strategies for 
the enhancement of thermal CP potentials developed in 
this work will thus present a valuable basis when consid- 
ering more complicated cavity geometries. 
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Appendix A: Scaling of potential depth with 
resonance order 

For a cavity of width a = vXko /2 the potential has v 
peaks, roughly located at z = — ^Afco/4 + (fi — l/2)Afco/2 
(p = 1 . . . v) and v — 1 minima at z = — ^Aj, /4 + fj.Xko/2 
(fi = 1 . . . v — 1). For a given resonance order v, the 
deepest minima are the ones closest to the cavity walls 
(and which have a maximum on both sides) , for example 
the rightmost one at z — {v — 2)Afco/4. It has to be 
compared with the lower of the two adjacent maxima, i.e., 



the one immediately to the right at z = (v — 3)A/co/4. The 
required depth of the deepest minimum is hence given by 

AU V = U[{v - 3)A fc0 /4] - U[(v - 2)A fc0 /4]. (Al) 

We calculate this depth for a cavity whose reflection co- 
efficients are independent of the transverse wave number 
k±, r p — — r s = r, and close to unity, i5 = l-r<l. Not- 
ing that the oscillating part of the potential is determined 
by C/pr and introducing the definition (|2"Tj) . cf. Sec. HID) 
we thus have 



AU v <x 1(1 



1v , 



(A2) 



We consider in the following only the terms which do 
not vanish as <5 — » 0+. For arbitrary <j>, we have 



2^3 A 3 fc 



Im 



2xu 2^/2 cos ^ 

ax — k—. (A3) 

1 - r 2 e lx v ' 



which, after expanding the fraction in powers of r 2 , solv- 
ing the integral over x and taking the imaginary part, 
can be written as 



Ofc j =0 



■ <f>) + yd 



with 



(A4) 



y(p) = - -3 + 



V 

— tt- sm(27wp) 

•p 1 



4Z/ 2 7T 2 
P 



cos(27rz/p) 



(A5) 



In particular, this implies 
2r 



H- - 



Ofc j =0 



i + i 



2/< 



(i + i-£) 3 (J + 1„? 



7(i-I)=-— V 



Ofc 



3=0 



r 2 :i 



3 + 1 



(A6a) 
(A6b) 



The evaluation of sums with simple denominators can 
be performed by using the relation (formula 9.559 in 
Ref. (H) 



00 r 2j 



E 

J=0 



1 



, = =-F(l,b;l + b;r 2 ), 
3+o b 



(A7) 



valid for any b ^ 0,-1,-2,... Here, F(a,b;c;z) = 
2-Fi(a, b; c; z) is a hypergeometric function which in turn 
has the following expansion in powers of 8 = 1 — r (for- 
mula 15.3.10 in Ref. (H) 



1 



F(l,6;l + b-r 2 ) ~ - ln£ - ln2 - 7 - t/>(&) (A8) 
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FIG. 10: Graph of <p(i>)/v (solid, red line) and its asymptote 
(dashed line) as given by Eq. (|A13[1 . 

as S — » 0+, the correction terms being of order <51n<5. 
Here, ip(x) is the logarithmic derivative of the gamma 
function and ^>(1) = —7 where 7 = 0.577216 is Euler's 
constant. 

For the sums in Eq. (|A6j) with cubic denominators, one 
can set r = 1 with an error of order 5. The sums are then 
simply Hurwitz zeta functions £(3, b) 

00 1 

£(7T^ C(3 ' 6) - (A9) 

1=0 v ; 

We thus find 

I{\ ~ -l)= -~[lnS + <p(u)] + ... (A10) 

vA ok 



as 5 0+ with corrections being of order 5 In d and 

+ ^[c(3,i-£) + c(3,£)]- ( An ) 

Some numerical values of tp{v) are 

ip(2) = - 0.1134423724; (Al2a) 
<p(3) = - 0.4015949503; (Al2b) 
<p(4) = - 0.7384479470. (Al2c) 

We give a plot of ip(y)jv in Fig. [TU] For large v, one 
obviously has 

re - 0.4091613938j + 0.9431471806 + ... (A13) 

plus terms of order and smaller. Already at v = 4 
this is a fairly good approximation to Eq. (|A11[) as shown 
in Fig. [10] 
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